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Abstract 

We derive the evolution equations for the spectra of a space ("Uni- 
verse"). Here the spectra mean the eigenvalues of the Laplacian on a 
space, which contain the geometrical information on the space. 

These spectral evolution equations are expected to be useful to analyze 
the time evolution of the geometrical structures of the Universe. In par- 
ticular, it is indispensable to investigate the time evolution of the spectral 
distance between two spaces, which is a measure of closeness between two 
geometries defined in terms of the spectra. 

As an application, we investigate the time evolution of the spectral dis- 
tance between two Universes that are very close to each other; it is the first 
necessary step for analyzing the time evolution of the geometrical discrep- 
ancies between the real Universe and its model. We find out a universal 
formula for the spectral distance between two very close Universes, which 
turns out to be independent of the detailed form of the spectral distance 
nor the gravity theory. Then we investigate its time evolution with the 
help of the spectral evolution equations. We also formulate the criteria for 
a good cosmological model in terms of the spectral distance. 



1 Introduction 

The notion of closeness or distance plays an essential role in physics. We fre- 
quently encounter the concept in the course of constructing and applying a the- 
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ory. 

The first step of finding out a law of nature is often classifying the objects 
in question into several categories (e.g., Hubble's classification of galaxies). For 
classification, we are tacitly assuming the concept of closeness between the ob- 
jects to be classified. Next, when we judge the validity of a theory by experiment, 
we construct a suitable parameter space associated to the theory and check the 
closeness between the point predicted by the theory and the experimental data- 
point in the parameter space. Here some kind of distance in the parameter space 
should be assumed. Finally, when we try to explain an observed phenomenon by 
a certain model based on a more or less established theory, we need to compare 
the observational data with the predictions of the model (e.g. the relation be- 
tween a signal of the gravitational wave and its templates; a comparison of the 
observational data with the results of numerical simulations based on a model). 
Again some notion of distance is required in a parameter space suitable for this 
purpose. 

Here we find out a universal setting for comparing "theory and reality" : There 
should be a parameter space equipped with the notion of closeness. We prepare a 
set of models (or templates). Each model corresponds to a point in the parameter 
space, so that the points corresponding to these models are distributed over the 
parameter space. Now, the observed data define another point in the parameter 
space. Then we try to find out the model-point that is closest to the data-point. 

Thus we realize the significance of establishing a parameter space and a dis- 
tance/closeness on it according to the problem we need to study. 

The same situation occurs in cosmology and spacetime physics. Indeed, we 
often need to consider a set of Universes, rather than just only "our Universe": 
Cosmology itself is a trial for grasping overall, averaged nature of the complicated 
reality of our Universe in terms of models; when we question why our Universe 
emerged rather than other possibilities, we are considering a set of Universes. 
According to the above considerations, thus, it is essential to establish a space of 
all Universes and a distance/closeness between any two Universes among them. 

In a series of investigations, it turned out that we can in fact construct a space 
of all compact Universes equipped with a sort of distance. 

Let Riem be the space of all {D — l)-dimensional, compact Riemannian ge- 
ometries without boundaries^. On Riem, we can introduce a measure of closeness 
in terms of the spectra, a set of eigenvalues of a certain elliptic operator [Q. Here 
we consider only the Laplacian A as an elliptic operator. For a given geometry 
Q G Riem, we get the spectra, or a set of eigenvalues of the Laplacian {A„}^q = 
{0 = Ao < Ai < A2 < ■ ■ ■ < A„ < ■ ■ ■ I 00}, numbered in an increasing order. 
Since A„ has dimension [Length"^], the higher (lower) spectrum in general re- 
flects the smaller (larger) scale properties of the geometry. Therefore the spectra 

^ Throughout this paper, D represents the spacetime dimension, so that the spatial section 
is always {D — l)-dimensional. 
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are desirable quantities for describing the effective geometrical structures of the 
space at each observational scale, e.g. the scale-dependent topology 0, Let us 
call this type of representation of geometry in terms of the spectra the spectral 
representation of geometrical structures 

Suppose Q and Q' are two spaces in Riem, and let {Xn}'^=o and {A'„}^o be 
the spectra for Q and Q', respectively. By comparing {A^j^^^ and {A'„}^^^ in 
a suitable manner, we can introduce a measure of closeness between Q and Q' 
of order IQ. It compares two geometries up to the scale of order 0(A^^^^), 
neglecting the smaller scale differences (we shall discuss more in detail on this 
topic inU and ^). 

However, it turns out that this measure of closeness dN^G, Q') does not satisfy 
the triangle inequality though it satisfies the other two axioms of distance 
Even though the triangle inequality is far from a must from the viewpoint of the 
general theory of point set topology, it is certain that the inequality makes several 
arguments concise and it makes the measure of closeness more compatible with 
our notion of closeness. 

This problem has been resolved by realizing that the breakdown of the trian- 
gle inequality is a mild one in a certain sense, and proving that Riem equipped 
with d^iQ^Q') forms a metrizable space 0. In other words, it has been justi- 
fied to regard dN{G,Q') as a distance provided that a care is taken when the 
triangle inequality is required in the argument; we also found out the alterna- 
tive of dN^G, Q') to be used when the triangle inequality is needed. (See and 
Ref. ^ for more details.) As an immediate consequence, we can even introduce a 
distance on Slice{J^, g), a space of all possible time-slices of a given spacetime 
{J^,g). (See Ref. |0] for more details on this point and its apphcation to the av- 
eraging problem in cosmology ^.) From now on, we call d^iG, Q') the spectral 
distance. 

Following the arguments at the beginning of this section, we have established 
a parameter space and a distance on it appropriate for spacetime physics: The 
space of all spaces of order A^, iSat, which is a completion of {Riem, d]y)/^ is what 
we had been searching for. (Here indicates the identification of isospectral 
manifolds . We discuss a physical interpretation of the isospectral manifolds in 
Section |. See also Ref. [|, |[.) Because of its nice property, Sn can be regarded 
as a basic arena for the study of spacetime physics. For instance, we can define 
integral over S^- [^], which would be essential in quantum cosmology. 

Being at hand, we are now in a position to handle a set of Universes. For 
definiteness, let us focus on cosmological problems now. In cosmology, we need 
to judge to what extent a model reflects the real Universe. There is no guarantee 
whether cosmology is possible, viz. whether a model close to the reality at some 
instant of time, remains so all the time. From the viewpoint of the spectral 
representation, this fundamental problem can be visualized as follows: Let Q be 
the real Universe at present with respect to a certain time-slicing. Let Q' be a 
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model located in the neighborhood of ^ in iSa? Then we should investigate 
the time evolution of diy{Q,Q') and should analyze in what conditions d]^{Q,Q') 
remains small during a certain period of time. 

This kind of investigation is now possible since the spectral distance is defined 
explicitly in terms of the spectra, which have a firm basis both physically and 
mathematically. What we now need is, thus, to analyze the time evolution of 
the spectra. In the spectral representation, the spectra {A„} are placed in the 
most fundamental position. Thus, from the purely theoretical viewpoint, too, it 
is interesting to investigate the time evolution of the spectra in detail. 

As a first step, we can understand how the time evolution of the spectra is 
induced by the evolution of geometry as follows: By evolving an initial {D — 1)- 
dimensional geometry (S, h) according to the Einstein equations (in the Hamil- 
tonian form if necessary), we get a 1-parameter family of geometries (S, 
In principle, then, we can get the spectra for each geometry (E,/i(t)). In this 
manner, we get a 1-parameter family of sets of spectra {A„(t)}. 

It is more preferable both theoretically and practically, however, if the time 
evolution of {A„(t)} is described (1) solely in terms of spectral quantities, without 
any explicit reference to the metrical information behind them, and (2) in the 
form of differential equations of {A„(t)} with respect to time. 

The key procedure for achieving this goal is to investigate the response of 
{A„} to the change of the spatial metric h. Since the latter is controlled by the 
Einstein equations, we thus expect to obtain the spectral version of the Einstein 
equations. The main aim of this paper is to obtain the fundamental evolution 
equations for the spectra, by putting this program into practice. 

In section ^ we prepare several formulas that are needed in the subsequent 
investigations. In section ^, which is the main part of this paper, we derive the 
spectral evolution equations. As basic applications of the results we obtained, 
we discuss three topics in section ||. In we study the spectral evolution of 
the closed Friedmann- Robert son- Walker Universe. In we study the spectral 
distance between two Universes that are very close to each other in iSat. We 
find out its universal expression in the leading order, which is independent of 
the detailed form of the spectral distance nor the gravity theory. In |4.3| , we 
investigate the time evolution of the spectral distance between two very close 
Universes. Section |^ is devoted for discussions. 

2 Basic formulas for the spectra 

Let (S,5f) be a (D — l)-dimensional compact Riemannian manifold without 
boundaries. We set an eigenvalue problem for the Laplacian A on {M.^g)^ 

^Here a model means a spacetime (usually it possesses much simpler geometrical structures 
than reality) along with a certain fixed time-slicing. We here regard an identical spacetime 
with different time-slicings as two different models. 
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A/ = -A/. Let {A„}„=o,i,2,.. := {0 = Aq < Ai < A2 < ■ • • < A„< • • ■ T 00} 
be the set of eigenvalues, or spectra hereafter, arranged in an increasing order. 
For simphcity of formulas, we assume that there is no degeneracy in the spectra 
throughout this paper. Let {/n}n=o,i,2,-- be the set of real-valued eigenfunctions 
that are normalized as 



(/m; fri) ■ / fm fn ^ ^mn i (1) 



where the natural integral measure on (S, g) is implied by y := y det((7a6). Let us 
note that a set {-^^t^a/nj^i forms an orthonormal subset of 1-forms on (S, (7), 

I dafm 9""' ^=Smn • (2) 

2.1 Spectral components of functions and tensors 

Let A and Aab be any function and any symmetric tensor field, respectively, on 
(S,(yf). It is useful to introduce diffeomorphism invariant quantities (A)^^ and 
{^ab)mn defined as 

(^ab)mn •= / ^ / ^"/m "^ab 9^ fn ^ ■ 

Here, for a quantity of type (^afe)m,n' always understand that n, m > 1 unless 
otherwise stated. We note that Eqs.(|l|) and (H) can be expressed as 

We also employ an abbreviated notation 

For later uses, we develop these notations to 

(^)lmn- ~ ^^fl)mn ' klmn ■~ i^f^) Imn ^ ' ' ' ' 

(Aabfk) 

Lmn ' 

We note that, for arbitrary functions A and B 

oo 

(AB)^^ = Y.{A)^,{B),^ . (4) 

fc=0 

To show this formula, we insert the 5-function into the integral expression of 
{AB)^^, noting that 5{x,y)r^ = T.T=o Ux) Uv) ■ 



The following formula is also useful, which relates {Agab) mn with {A) 

{Agab)mn = r, /T T={^)mn + r, / a a i^^) mn ■ i^) 

To show this formula, one modifies the defining equation for {Agab)^^ with the 
help of partial integrals, noting that the R.H.S. (right-hand side) should be sym- 
metric in m and n . 

Setting m = n in Eq.(|), we get 

{Ag^b)^={A)^ + ^JAA)^ . (6) 
2.2 The variation formulas 

We frequently consider the variation 6Q of a certain quantity Q below. Here we 
treat S as a general variation for the time being. In later applications, the time- 
derivative (the Lie derivative along a time-flow vector in a spacetime picture) 
is mostly considered as the variation operator 6. 

Now, noting that A/ = ^da{^ g"'^df,f) for an arbitrary function /, the vari- 
ation of A is represented as| 

6Af = ^da{g ■ 5g)d^f - -da{g''%bcdj ^) . 
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Thus, employing the same kind of notation as (^4)^^, we can introduce the 
quantity 

(5A)_ := I U6Af^^=y Udaig ■ 6g)d^fn^- J fmdaig^'SgbcdJn^) , 

where we note that the variation is taken only for the operator A, and not for the 
eigenfunctions /„, fm- We should also keep in mind that {6A)^^ is noi symmetric 
in m and n, unlike {A)^^, because A is an operator, and not a function. Noting 
that /o is a constant function (see Appendix^, it is evident that {6A)^q = 
and (SA)^^ = ^{g ■ 5g),^ (m = 0, 1, 2, ■ • ■)• 

With these preliminaries, we now investigate the variations of spectral quan- 
tities. We start with the variation of the spectra. From Eq.( |BlCI| ) in Appendix^ 
it follows that 

5K = -{6A)^ , (7) 
which is a basic result of the perturbation theory ("Fermi's golden rule"). 

^ To avoid trivial indices, we flexibly adopt notations such as u :— u"-, A ■ B :— AabB'^^, 
u ■ V := Ua and A.' := A^ . We also flexibly choose symbols for the derivative of a function 
/, such as Daf = Df = daf- (Here Da and D denote the covariant derivative.) 
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Now let us investigate the variation of the eigenfunctions. We have a gen- 
eral formula Eg . ( P 1 1|) with Eg. ( P12D (see Appendix^ for the perturbation of 
eigenvectors. Here we need to specify the factor c^^-* in Eg. (|B12| ). For this pur- 
pose, we take the variation of the both sides of Eg.(|l]) for m = n. Noting that 

= 1^ ■ ^9^y, we get 

ifn, SU) = -^{g.6g)^ . (8) 

In the standard perturbation analysis in guantum mechanics, the inner-product 
is fixed, and not perturbed, while the eigenfunctions are perturbed. Thus, Sfn 
should be perpendicular to if is normalized. In our case, on the other hand, 
the inner-product is also subject to the variation because of the presence of the 
integral measure y. Thus 6fn is not in general perpendicular to /„, as is clear 
in Eg.(D. Combining Eg.® with Egs.fpTlD and Eg.([BT2D, we see that in 
Eg.( |B12| ) should be chosen as c^^ = —\{g- Sg)n the present case. Thus, we 
get 

oo 

^fn = J2fk f^kn , (9) 
k=0 



where 



for m ^ n 



I -l{g ■ og)m for m = n . 
Taking the inner-product of the both sides of Eg.® with /m, we get 

/imn = (/m, , (H) 

which gives a clear interpretation of the guantity Hmn as the projection of 5/„ to 
the direction of fm- 

With the help of Eg.®, it is now straightforward to show that 

^{^)mn = {S^)mn + ^ " ^^?)™„ + T.i^) mk l^kn + E(^)nfc l^km • (12) 

^ k k 

In particular, for the case of m = n, we get 

HA)^ = {5A)^ + ^{Ag.5g)^ + 2j:{A)^,l^,^ . (13) 

k 

Introducing Tmn ■= -l^mn -\{9 ■ Sg).mn (^^^^ ^^at Tnn = 0), Eg.(0) is also 
represented as 

k k 

It is interesting that this expression is in a similar form as the covariant derivative 
of a symmetric tensor. 
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In the same manner, we get a formula for ^{Aab)^^ as 



'5(^afe)„„ = (SAab),^^ - (5 In ^A„A„ ){Aab)^^ + -{Aab 9 ■ Sg)^^ - {A^'Sga + A^%g,a)mn 

k ' k ' ™' 

and for the case of m = n, we get 

HAab)n = {SAab)r^-{SlnX^){Aab)^ + ^{Aabg■6g)^-2{{A6g)ab)r. 



+2j2J-^{Aab)^,fikn . (15) 
k V 

2.3 Basic identities 

We obtain important relations by taking the variation of the basic identities 
Eq.(0). 

First, we take the variation of the both sides of 6mn = {^)mn- Then, with the 
help of Eq. (|12D , we get 

which implies an identity, 

■ S9)nin = -2(/imn + fJ-nm) ■ (16) 

Thus, with the help of Eq. (p!oD , we get 

{9 ■ 89) = -Y^^^^^^rnn ' W nJ ^OT TU ^ U . (17) 

On the other hand, no condition is imposed on {g ■ Sg)^ except for 

{9-Sgh = 2^-^ , (18) 

which follows by an independent argument (see Appendix ^ . 

Now we take the variation of the both sides of 6mn = {9ab)mn m,n > 1. 
Then, with the help of Eq.(|^), we get 

= -(<5fi'a6)m„ - ('^lnA„)5m„ + Wy^ + W /i„^ , (19) 

V V '^m 

where Aab '■= Aab — \A.'gab for any symmetric tensor Aab- 
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A m ' 

- Am-^ (v^ - y£ (5A)_) for m ^ n 

— 51nA„, — \{g-^g)n for m = n . 

(20) 

The fact that (5A)^q = (m = 0, 1, 2, ■ ■ ■) suggests us to formally define that 
{^9ab)mo = (m = 0, 1, 2, ■ ■ ■). We adopt this formal definition here for a nota- 
tional neatness (see Eg . (^1]) below). 

From Eqs. (p!7D and (pO|), we obtain a formula for {6A)^^, 



We can also derive Eg . (^I]) directly from the definition of (5A)^,^. 

We note that Eq.(^) is valid for the case of m = n also, due to Eq.(|^) and 
the second equation of Eq.(^). Furthermore, we realize that this formula is 
also valid for m = or = on account of the formal definition {5gab)mo ~ ^ 
(m = 0,l,2,---). 

We now investigate a different type of identities. 

We pay attention to the second equation in Eq. (pOD , 



5\nK = -{5gJ^^-\{g-5g)^ . (22) 
From this relation, it is straightforward to show a formula 

^a.^ = -(A.,)„-i(A-)„ , (23) 
(>gab ^ 

where Aab is any symmetric tensor field. 

Let u"" be any vector field. Substituting Aah = gab, the L.H.S. (left-hand 
side) of Eq. (p3| ) vanishes because of the diffeomorphism invariance of the spectra, 
so that 



{Cu gab)n = {D-U gab )„ - " «), 



(DiaUb))^ + -{D-u)^ = . (24) 

We now note another identity. Using the basic properties of the covariant 
derivative Da, it is easily shown that 

ADaf = DaAf + RabD'f , 

where / is any smooth function. In particular, choosing /„ as /, we get 

igabA-Rab)d'fn = -Xndafn , (25) 
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i.e. dafn turns out to be an eigenfunction of Qab^ — Rab with the eigenvalue A^. 
Taking the inner-product of Eq.(^) with dafmi we get 

I DaD,fmD''D'U^ = -^/KJ'n{Rab)^^ + Xl5^n , (26) 

or 

{Rab)^n = {9ab^)mn + ^n^mn ■ (27) 

3 Evolution equations for the spectral quanti- 
ties 

Now we investigate the time evolution of the spectra of a space. We have prepared 
in the previous section basic formulas regarding the responses of the spectral 
quantities with respect to a change 6gab- When we let 6gab be of a dynamical 
origin, thus, we automatically obtain the evolution equations for the spectral 
quantities. 

We consider (E, h), a (D— l)-dimensional compact Riemannian manifold with- 
out boundaries, as a mathematical model of the spatial section of the Universe. 
For the present purpose, we interprete a quantity 5Q as ^| ^^Qia), and identify 

the latter quantity with Q := C^; Q, where t is a time-flow vector]^. In particu- 
lar, we replace Sgab, g ■ Sg and 5(7^^ in the previous section with corresponding 
quantities as 

Sgab I > hab = 2NKab + 2D^aNb) , 

g-6g I — > h^^hab = 2NK + 2D-N , 
Jgab ^ Kb = 2NKab + 2'D^) . (28) 

Here Kab is the extrinsic curvature and K := K, '; N and Na are the lapse function 
and the shift vector, respectively. 

3.1 Evolution equation for {A^} 
First, Eq . (^2]) becomes 

A. = -{2{NK,,)^ + {NK)^)X^ 

= -2(iVir,,)„A„ + i(A(iVir))„ , (29) 

where we used the identity Eq.(|2^) in the first line, and Eq.(^ in the second 
line. We note that the shift vector A^^ does not appear in Eq . (pOf) , reflecting the 
spatial diffeomorphism invariance of A„. 

^ See the argument at the end of Appendix H also. 
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Now, Eqs. (p!6|) and ( pOD becomes 
{NK)^^ + {D ■ iV)^„ = + 

where /i^^ is given by Eq.(|OD with ((5A).^„ replaced by (A)^^. 
Eq.(pTD becomes 

(30) 

Ultimately we are only interested in the dynamics of the spectra {A„}, which 
is invariant under the spatial diffeomorphism. The coordinate dependence of 
other subsidiary variables like ^mni which can dependent on the shift vector Na, 
should not have any influence on the dynamics of {A„}. Hence we can set A^a = 
from the very beginning to make formulas simpler. For simplicity we also set 
N = 1 hereafter. 

Now let us investigate the evolution equations for the spectra, Eq . (p9D , in 
detail. At this stage, it is useful to develop notations. First, we note that we can 
expand any scalar function A{-) in terms of {/n}^o- 

oo 

A{-)=Y.^nU-) ■ (31) 
1=0 

Then, it follows that 

oo 

AAi-) = -Y^KAr,fn{-) , 

1=1 

where we see that the homogeneous component of A, Aq, does not appear in the 
summation. 

The homogeneous component Aq is related to the spatial average of A over 
the spatial section S, A^v '■= y A^, as 

A,^ = Ao/Vv . (32) 
Next, let us introduce the quantity 

(/ m n) := = = J fi fm fn^ ■ (33) 

Note that {I m n) is totally symmetric in /, m and n. We can also introduce 
a similar quantity {I m n k) := J fi fm fn fk y/ , and similarly {I m n k h), 
and so on. However, the quantities of the type (/ m n) are sufficient since other 
quantities can be expressed in terms of {I m n). For instance, 

{I mn k) = -J2{{l m n k) + {m n k I) + {n k I m)} , 
3 v 
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due to Eq.(|). 

We note that the quantity of the form (^) is represented in terms of Ai 
and (/ m n) since 

{A),mn- = EMl'lmn ■■■) . 
I' 

In the same manner, the quantities of the form (Aab) ki... represented in 
terms of {Aat) i ^„ and (/ m n) because of the relation 

V 

In particular, the quantity of the form {Aab) can be represented in terms of 

Keeping the applications to cosmology in mind, it is also useful to introduce 
the quantities eab and Tab defined as 

Cab := Kab - J-) _ ^ Khab , 

Tab ■■= Rab- \ ^ Rhab ■ (34) 

Note that e. ' = r. ' = 0. The quantities eab and Tab characterize the deviation of 
the spatial geometry (S, h) from the isotropic geometry. 
Now, we note that 

and 

D — 3 

{Kab)mn ~ ~ 2[D — 1) ~'~ ^^"*^mn 

Here Eq.(^l]) has been applied to K and we have used Eq.(|^). 
Thus, Eq. (p9D is represented as 



A 



2 



1 E(An + -^h) Ki {lnn)~ 2Xn{eab)^ • (35) 



This equation is also valid for n = 0, viz. it is compatible with Aq = 0. 

Looking at the R.H.S. of Eq. (p5D , it turns out that we further need the equa- 
tions for (/ m n) , Ki and {eab)n ■ We investigate them one by one below. 
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3.2 Evolution equation for {I m n) 

Applying Eq.(|l2D along with Eq.® to = {I m n),we get 

[l mn) {(/ m I') fii>n + [jn n I') pm + {n I I') fii>m} + ^{l m n l')Ki> . 

V v 

(36) 

Here, from Eq.(|10]) with Eq.(|g) (A^ = 1, A^^ = 0), pmn is given by 

""^^ ' E.(An-fef(A„-AO)KK/mn) + ^^(e,,)_ for m 



2{D~1) A„-A„ 



mn 



Klil m m) for m 

(37) 

3.3 Evolution equation for Ki 

Now, since Ki = {K, fi), we get Ki = {K, fi) + (AT, fi) + (AT, fi AT), resulting in 
it; = {k)i + ^ + E ^" , (38) 

I' VI" 

where Eq.(|lT]) has been applied. To get the detailed expression for the first term 
(A');, we recall basic formulas of the canonical Einstein equations. We note 

/<■ = -i^p - ^^AR + A- + . A-) + . (39) 

2a D- 2^ 2(A) -2)^ A) - 3 ' D-2 ' ^ ' 

along with the constraints 

i?-A:-A: + A:2-ip-2A = , 

a 

D'Kat-DaK+^Ja = , (40) 

2a 

3 

where a := j^;^ and A is the cosmological constant. Here we define the quantities 
p, p, Ja and Sab in connection with the above equations: With the help of the 
normal unit vector of the spatial section, the energy-momentum tensor of 
matter T"^ can be decomposed into three components: = +T^^'^^ +T^^^^ , 
where each term has the following form, 



XX 



(The suffix A implies "perpendicular to the space" and || imphes "along the 
space".) Here J° and S"-^ are spatial quantities and they are uniquely identified 
to their spatial counterparts, J" and S'"'', respectively. Then p, J a and Sab are 
interpreted as the energy density, the momentum density and the stress tensor of 
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matter, respectively, and p := j^S^"" defines the pressure of matter. (Tlie vector 
and tensor indices of spatial quantities in this context are lowered and raised by, 
respectively, the spatial metric hab and its inverse /i"^.) 

Now, Eq.(p9|) can be modified by means of the first constraint equation 
(Hamiltonian constraint) in Eq. (^OD . In particular, the following two forms would 
be useful for our purposes. First, by eliminating R from Eq.(^), we get 

K = - — - — -(p + - -p)-K-K+^^A . (41) 

Second, by eliminating the term K ■ K from Eq.(^), we get 

A- = ±^(p-p)-A--fi+?(^A . (42) 

Thus, we obtain the equation for {K)i based on Eq. (^T]) , 

(43) 

In the same manner, based on Eq.(E2[), we get 



We also note that, from Eqs.(|T8|), (|3|) and ( [X2D , 

'^ = K^|^/V = K^, . (45) 
The first constraint equation in Eq.(^) is translated into 

+ E KvKv\l' I" I) - -Pi - 2AVv 6oi -{e-e)i = . (46) 

U - I'pjf, a 

We also note that, taking the inner-product with j-^Dafi (/ 7^ 0), the Bianchi 
identity D^Rab — \DaR = turns to 

J2^ + ^^^^P'^^V,,), = (1^0) . (47) 

Now, taking the inner-product with j-^afi (^ 7^ 0), the second constraint 
equation in Eq.(^) (momentum constraint) turns to 

+ + A}=0 (/^O) . (48) 
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Finally, noting that 



cd, 



we get 

= ^ ^" ^) + E Ri'f'i'i --{D- J)i - 2(e ■ r), , (49) 

where we used the momentum constraint Eq.(^) to reach the final form. 

3.4 Evolution equation for {eab)i^rnn {rab)i^mn 
First, we derive the evolution equations for Eab and Tab {N = 1, Na = 0), 

1 D — 3 

eab = —{Sab - phab) " 7:; — 'Keab - Tab + 2(e ■ e)ab , (50) 
2a U — 1 

Tab = --^{AKKb + iD-3)DaDbK} + --^D-Jhab 
D — 1 a D — 1 

2 2 

+ 7;^ -e ■ rhab ■ (51) 

From Eq. ([50|) along with Eqs.(^ and (p^ , we get 



k ' fc ' ™ 

,mn 

k 

In the same manner, we get from Eq. (pT]) , 

("^abii^rnn = " (in \/A^A„] (r„b); - ^ ^ Ki,{l' I I") {Tab) I. 

"'"X] \/ '\~i''^ab) i^rnk^^kn + \/ T (''^a.b) i^nk^^km 
k ' fe ' ™ 

,77171 ' 

(53) 

k 

where 

1 D — 3 



I'l" 
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1 1 



a 



^ J- iiiii ^ J- inii 



^ ^ I'l" 



(54) 



For convenience, we also present the formulas Eq.(p2D and Eq.(^) especially 
for / = 0: 



1 



^ V y!j V '^m 

-2((e- e)ab)^„ • 



(55) 



(rab), 



where 

(rab), 



{rab)mn- In V A™A„ (r^b) 



_5-D 



+ X] V \ {''"ab) mkt^kn + X] V ^ {^ab)nkl^kr, 



-4(e%rb)e), 



1 — 3 

^ _ YXkKk{hab)k,mn " _l Y ^k{DaDbfk) mn 

k k 

1 1 - - 2 

^ ^ ■ J)k{hab) k,mn " ^ _ ^ ^ki^ab) k,; 



aD 



(56) 



-(Ae„,)_ + 2(Z}^D(,e,),)_ + ^—j ^^(e ■ r)fc(/i,,),__ . (57) 



Eqs. (|5D-(|8D, (H) (or (0)), (|5D-(||), and (|5|)-(||) are the fundamental equa- 
tions for the investigation of the spectral evolution of the Universe. 

They form hierarchy equations and we can continue to get equations for higher 
hierarchies. This hierarchical property is a reasonable consequence since we are 
looking at the global quantities, and not the local ones. In practical applications, 
thus, we need to make a suitable truncation. Typically, we get equations of 
higher order in e^fe and Tab when we continue to go into further hierarchies. We 
can often regard eab, rab, and their spatial derivatives are small in cosmological 
applications. In such cases, the truncation becomes a reasonable approximation 
procedure. 
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4 Basic applications of the spectral equations 

4.1 The Friedman- Robertson- Walker Universe 

As a basic example, let us consider a closed Friedmann-Robertson- Walker Uni- 
verse. 

We set eab = and r^fe = 0. We also set J = 0. From Eqs. (g^) and (HD, we 
get Ri = and Ki = (1^0). Then, we get = (/ ^ 0) from Eq.®, noting 
that (0 0/) = ^5o«. Eq.(|Sp along with Eq.(||) imply that pi = (l ^ 0), 

noting that = —^5mnKo/VV in the present case (Eq.(P7D). 
Now, Eq.(^ with Eq.(^) yield 

X y 

thus, 

It is a simple scaling behavior expected from the dimensionality of A„. 

Noting that /i,o = -|5/o^o/v^ (Eq-d^)), we get from Eq.® for / = with 
Eq.®, 

Thus, 

Ro{t) = (^^j Ro{0) , 

or, noting Eq. (|32|) , 

It is also a simple scaling behavior expected from the dimensionality of R. 
On the other hand, from Eq.(^6D for / = with Eq. (^Sj) , we get 

D-2 fV\^ 1 

—l[v) +«.v--P„-2A = , 

or in a more familiar form, 

'd\2 k 1 1 



aj {D - 1){D - 2) a^'^ {D-l){D-2)^ ° ' 

where we introduced the scale factor a as ^ oc a^^^ and the curvature index 
k := i?av(0)a(0)2. 

It can be solved once the matter content is specified. 
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4.2 A universal formula for the spectral distance between 
two Universes that are geometrically close to each 
other 

The spectral evolution equations developed in the previous sections are of essen- 
tial significance for the study of spacetime physics along the line of the spectral 
representation of geometrical structures. 

We have the spectral distance dN{G,G') as a measure of closeness between 
two spatial geometries Q and ^' It measures the difference between Q and Q' 
as the difference of "sounds" of them, i.e. the difference of the spectra. 

Then we are given a basic arena for spacetime physics, i.e. the "space of all 
spaces" iSat; it is basically a space of all compact Riemannian geometries equipped 

with d^{g,g') i. 

Since we now have the spectral evolution equations developed in the previous 
sections, we can investigate the time evolution of dN{Q-,Q')- The fundamental 
importance of such an investigation becomes clear when we take Q as the real 
Universe and Q' as a model Universe which is expected to be "close" to ^ @! 01- 
There is no guarantee that the model Universe remains a good model for the 
real Universe in the future also P, ^. Stated in terms of the spectral distance, 
there is no guarantee that di\i{Q,Q') remains small in the future. The detailed 
investigations along this line would be done separately. Here we only look at 
some basic features of the evolution of dN{G,G') when Q and Q' are very close 
initially in iSat. 

In this subsection, we derive a universal formula for dN{G, G') when Q and Q' 
are geometrically close in Sn-, which is valid independently of the detailed form 
of dtqiG^Q') nor the gravity theory. In the next subsection, we discuss its time 
evolution. 

First, it is appropriate to recall basic settings of the spectral representation 
We set the eigenvalue problem on each manifold Q and Q' , 

^f = -\f , 

then the set of eigenvalues (numbered in increasing order) is obtained; {Xm\m=o 
for g and {A'J- o for g' . 

Now the spectral distance between g and g' is defined as @, 0] 

dM{g,g') = j:^{Y'] ' (58) 

n=l \^n/ 

where J-'{x) (x > 0) is a smooth function which satisfies > 0, J-'{l/x) = J-'{x), 
T{y) > J^{x) if y > X > 1 and J^(l) = 0. Then, it follows that J^"(l) > 0. 
However, in order to let d]^{g, g') detect a fine difference between g and g' when 
they are very close to each other in S^, we further postulate that JF"(1) > 0. 
(SeeEq. (^.) 
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In particular, it is convenient to choose as JF, jFi(x) = ^\n^{y/x + l/^/x). 
Then we get 

Ma.ao4|;.ni(/J+/|) . (59) 

It turns out that d^iQ^Q') does not satisfy the triangle inequality. However, 
it does not cause a serious problem. Let Riem be the space of all {D — 1)- 
dimensional, compact Riemannian geometries without boundaries. Then it is 
proved that the space [Riem, d^) /~, where d^ is given by Eq.(|59D, is a metrizable 
space 0. (Here ^ indicates identification of isospectral manifolds 0.) It justifies 
to regard d^iQ, Q') as a distance, provided that we are careful when the triangle 
inequality matters in the argument. 

The above property is shown in the following way. Since the breakdown of the 
triangle inequality turns out to be a mild one in a certain sense 0] , it is expected 
that one can make a slight modification of to recover the inequality. Indeed 
we can find ^^0(2^) •= |lnmax(A/x, as a modification of J^i. In this case, 

Eq. (|58|) becomes 

1 ^ 

dniQ, Q') = T^Yl In max 

^ n=l 

It is easy to show that dN{Q,Q') satisfies all of the axioms of distance, so that 
it is a distance. Now, one can show that {Riem,di\f)/^ and {Riem,di\i)/r^ are 
homeomorphic to each other. Thus, {Riem,dj\f)/^ is a metrizable space since 
{Riem, d^) / is a metric space. Let Sn be a completion of {Riem,dN)/^- The 
definition for dN has a more convenient form than the one for dN since the latter 
includes max symbol. Thus, in the actual applications, d^ is more useful than 
Jat. On the other hand, when we need to discuss mathematical properties of Sn 
precisely, d^ is appropriate. 

Let us consider the situation that Q and Q' possess the same topological 
structure and that they are very close in Sn- One can imagine that Q' = (E, h') 
represents the real Universe at some instant of time, and Q = (S, h) is a model 
Universe corresponding to Q' = (S, h'). We introduce the difference in the spatial 
metric 

7afe := Kb - hab , (60) 

and we treat jab as a small quantity. 

We regard the model Universe as a reference point, based on which we evaluate 
several quantities. In particular, we make use of the time-slicings of the model 
and investigate the time evolutions with respect to them. 

Now, from Eq.(|2^), we get 

S^^^n = ^^^^^ = -{lab)n-l{l)n , (61) 
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where 7 := h'^'^jab- 

Looking at Eg . (|58D , we see that 



An 



Here we note that = ^'(1) = 0, and e indicates a small guantity in the 

same order as 7. 

Leaving only the leading term, we thus get 



1_ A/, , 1 



n=l 



dN{g, G') = 2-^"(l) E [{lat)n + 2(7)n ) ■ (62) 



Here we note the postulation JF"(1) > 0. 

It would be also helpful to view Eg . (|62D as 



d^{g,g') = ^r'{i)i-i , (63) 

where 7 is a vector in R'^ whose n-th component is (7afe)„ + |(7)„; and a standard 
Euclidean inner-product is implied. 

Hence, we get a universal result on dN{g,G') when Q and Q' are very close 
in (Sat: The leading behavior of the spectral distance d^{Q, Q') is given by Eq. ^6^ ) 
( or Eq. j^UBj ) ), irrespective of the detailed form of the spectral distance nor of the 
gravity theory. 

In the case of Eg. (|59|) , we have chosen as JF, J^i{x) = ^\n^{^/x + l/^/x), 
hence JF"(1) = |. Thus, we get 

1 A /, , 1 



dNig,Q') = -j:[i^ab)n+2^l)n) , (64) 



n=l 

or 

dN{g,G') = ^7-7 , (65) 
where 7 is the same vector as in Eg. 

4.3 Time evolution of a small geometrical discrepancy be- 
tween the real and a model Universes 

Here we investigate the time evolution of the spectral distance for the two 'nearby' 
Universes as described in the previous subsection. Taking the time derivative of 
the both sides of Eg. (^) , we get 

MS,S') = -^E^Ainf\ . (66) 
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Now Eg. (p5|) can be represented as 



(In A„)- = --^ E(l + ^ ^) - 2(^a.). . 

In the cosmological problems, it is often useful to separate the term for / = from 
the terms for / > 1 in the summation like in the above equation: This separation 
is useful when the Universe is described by a homogeneous geometry plus small 
perturbations. Noting Eq.(E5|), we thus get 



(67) 



Let us define 



H 



V 



D~IV 

1 D-3A, 



-,T,(l^—j-)K,ilnn) 



^ l>l ^ 

Here H is an analogous quantity to the Hubble constant; in is attributed to the 
inhomogeneity while an is to the anisotropy of the geometry. Then, the formula 
can be represented more concisely as 



(lnA„)- = -2if„ , (69) 

where 

■=H + i„ + an . (70) 

The quantity Hn can be interpreted as the effective Hubble constant observed 
at the scale A~^/^ since it determines the rate of change of A„. Thus, Eq.([70|) 
describes the modification of the effective Hubble parameter at the scale A~^/^ 
due to inhomogeneity and anisotropy of the Universe at that scale. 

Leaving only the leading terms in Eq.(|66|) with the help of Eqs. (|69|) and (|70D , 
we thus get 

MQ.Q') = \i-5H , (71) 

where 5H denotes a vector in whose n-th component is 5Hn ■= H'^ — 

On the other hand, from Eq.(pPD with Eq. (pT|) , we can derive 7 = 26 H, which is 

compatible with Eqs.(|65|) and (ffll). 

In this subsection, SQ denotes the difference in a quantity Q for two spaces G and G' , 
defined as Q for G' (the second entry of div(-, •)) minus Q for G (the first entry of dN{-, •))■ 
Wherever necessary, we employ the notation 6{ ■ } (rather than 6{ ■ )) to avoid any confusion 
with the 5-function. 
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From Eg. ([7I|) , we get 



Mg,g') = ^6H-6H+^'y-6H . (72) 

Let us introduce g„ := —(1 + ■^), which can be interpreted as the effective 
deceleration parameter at the scale A~^/^. Then, Hn = —(1 + qn)H^, so that 

6Hn = -2{l + qn)HJHn-Hl6qn . (73) 

Since 5Hn appears in Eg. ([7^) only in the form of 7- it suffices to estimate 
SHn by leaving only the leading terms in Eg. (|73D . First, using Eg.([70|), we note 
that 

2 2 1 

gn - (1 - J^{tn + Oin))q " J^i^n + Oln) " ^(4 + "n) , 

where g := — (1 + j^). 

Next, it is also straightforward to get estimations 



^1 + qn)HnSHn 



[1 + q){H - {tn + Oln)} - -^(4 + dr. 



6H+{l + q)H6{Ln + 



Hl5qn ~ H^5q + 2{{l + q){Ln + an) + ^{Ln + dn)]5H 

ri 

-2{l ^ q)H 6 {in + an] - 5 {in + dn} ■ 

Looking at Eg. (|73D , we thus estimate 

5Hn--5{{l + q)H^} + 5{Ln + ani . 

Finally, we reach the estimation 

1^-1^ 1 
d^{g, g') -6H .6H--{Y: 7n)5{(l + q)H'} + -7 ■ 6{i:+ a} , (74) 

^ n=l ^ 

where 7„ := {^abln + |(7)n n-th. component of 7, and t and 6l are vectors 

whose n-th components are t„ and respectively. We note once again the 
definition of 5Q, the difference of Q in ^ and in this subsection (see the 
footnote just after Eg.([7T|)). The difference in metric, 7afe, is also defined in the 
same manner (Eg.(^)). Thus, Eq.(|7^ is symmetric in g and g\ as it should be. 

To make a detailed study, we need to investigate {eab) n ^'^d {vab) n ^Iso. It is 
helpful to note that /i^n in Eg.(P7|) gets simplified in the present case as 

1 V 

= -^^H6n,„, + 0ie) . (75) 
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Now, we set m = n in Eg. (|55|) . First, let us omit the first and the last terms 
on the R.H.S. of Eg . (|55|) . Next, we leave only the k = part in the third term, 

which gives -^jy{^ab)n = -'^H{eab)n ^i^^ ^^^^P Eg.(||). Next, with 
the help of Eg. (|75D , the fourth and the fifth terms on the R.H.S. in Eg.(^) are 
approximated together as —y{£ab)n = ~ ^)H{eab)n- Finally, with the help 
of Egs.(l^) and (|70D , the second term in Eg. (|55|) yields 2H{eab)n within the order 
of magnitude discussed now. Thus, we estimate 

(ea.L^-p-l)^M„-(ra6)„ . (76) 



In a similar manner, we can estimate {r'ab) n t'Y Egs. (|56D and (|57|). We set m = n 
in these eguations. Among all the terms in Eg.(^), it is a reasonable approx- 
imation to leave only the = part of the fourth term on the R.H.S., which 
gives —j^^{^ab)n- (This approximation is especially effective when we dis- 
cuss the large-scale behaviors (low-lying spectra) such as A~^/^ >> c^r, where Cg 
is the sound velocity of the matter and r is the typical cosmological time-scale 
of interest.) In Eg.(|56|), we omit the last term on the R.H.S. The third term 
can be estimated as —^^^(Tab)^ = ~(5 — D)H{rab) while the fourth and the 
fifth terms together are estimated as — ^(?"afe)„ = —{D — l)H(rab)n- Finally, the 
second term yields 2H{rab)n- Thus, we can estimate as 

2 

{rab)n - -2H{rab)r, " JJ^R^^A^ab) n ■ (77) 

Now, we can estimate dN{G,G') further based on Eg. ([7iD . First, d„ is given 
by Eg.(|7|), 

--(/}- l)i7«„-(r,fe)„ . (78) 

Second, regarding t'n, we go back to the definition of in (Eg. (|68D ). With the help 
of Egs. (|36D , (|38|), (|43|), (|69|) and (^), it is straightforward to get estimations 

0{e) , 

^~^H{lnn) + 0{e) , 




2 



-2Hl„-^ — ^-Mn , (79) 



where 



1 r/ \ . ^ D-3 1 , . ^ 
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With the help of Eqs. (fTHD and (|7^), Eq.([71D can be expressed as 



1 - - 1 ^ 
d^{g, Q') ^ -5H ■5H--iY: 7n)5{(l + q)H'} 

-\'f-s{^HZ+iD-l)Hd + c+^^:^^^M^ , (80) 

— * 

where c and Ai are vectors whose n-th components are, respectively, {rab)n 

Mn- 

Eqs. ( |7TD and ( PDD give us an idea on the factors that influence the 

validity of a cosmological model. 

First, a particular combination 7„ := {jab)n + determines the spectral 

distance d]si{Q,Q'). Second, 7 ■ 5H determines the rate of change of dN{G,G'), 
dN{G,G')- Whether d]<q{Q,Q') decreases or not is governed by the relative direc- 
tions of 7 and 5H in R^: One of the criteria for a good cosmological model 
would be that it makes the quantity 7 ■ 5H negative, or non-negative and small 
at least. Third, the acceleration, d]y{Q,Q') is determined by several factors: The 
difference in the effective Hubble parameter, 6H, has always a repulsive effect. 
Even though 7 = initially, so that (i^v = rf^r = initially, the spectral distance 



increases if 6H 7^ 0. It is required that the other terms in Eq. (|80D as a whole 
should be negative or at least non-negative and small in order to get a good 
model. 

Let r be the typical time-scale with which we want to discuss the evolution of 
the Universe. Then we can summarize the criteria for a good cosmological model 
as follows: 

(CI) dN{Q,Q') is small. 

(C2) TdN{G,G') is negative, or at least, non-negative and small. 
(C3) r'^d^iQ is negative, or at least, non-negative and small. 



5 Discussion 

In this paper, we have derived the spectral evolution equations of the Universe. 
A set of these equations forms one of the essential elements of the general scheme 
of spectral representation |jl[]: Now we have a space of all spaces 0, iSat, equipped 
with the spectral distance c/at, and the evolution equations on Sn- 

The spectral evolution equations are expected to be useful for studying the 
time evolution of the global geometrical structures of the Universe, since the 
spectra are especially suitable for describing global properties of a space. 

The significance of the spectral evolution equations becomes prominent in 
situations when we need to handle a set of spaces rather than just one space, e.g. 
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the comparison between the real Universe and its model, the relation between the 
underlying topological structures and its low energy behavior (scale-dependent 
topology 0, 0), and so on. 

As one of the important applications of the spectral evolution equations, we 
can now investigate a fundamental problem in cosmology: Whether cosmology is 
possible, viz. under what conditions and to what extent a model reflects the real 
spacetime faithfully ("The averaging/model-fitting problem in cosmology P, ^). 

As the first step in this direction, we have investigated in §^ the spectral dis- 
tance between two very close Universes and its time evolution. It is interesting 
that the spectral distance in this situation is universally determined by the quan- 
tity 7„ := (7afe)„+|(7)„ each scale n, irrespective of the detailed form of the 
spectral distance nor the gravity theory (Eq. (|63|)). Even though it is a special 
case when two geometries are very close to each other in Sn, this result would 
still provide us with a basic understanding about what kind of geometrical dis- 
crepancies would matter in the context of the validity of a cosmological model for 
the real Universe. More systematic studies along this line are certainly required, 
which would be presented elsewhere [p!0|] . 

Finally, we make some remarks on the spectral scheme in general. It is a 
different way of viewing geometrical structures from the standard way of de- 
scribing them. In this scheme, the geometrical information on a space is rep- 
resented by a collection of the whole of the spectral information measured by 
all available elliptic operators on the space H, 0. In ordinary cosmological 
observations, we use a particular observational apparatus so that we naturally 
get only a portion of the whole geometrical information on the space, according 
to which apparatus (or mathematically, which elliptic operators corresponding to 
the apparatus) has been utilized. Thus, it can happen that the geometry of the 
space cannot be fully identified by using just a single apparatus (some particular 
elliptic operators). This consideration provides us with a physical interpretation 
of the isospectral manifolds. From this viewpoint, there is no surprise in the 
existence of the isospectral manifolds. The spectral scheme describes the scale- 
and apparatus-dependent geometry of a space in a natural manner. According 
to what this scheme suggests us, there is no absolute model for the real Universe, 
rather a good model for reality depends on the observational scale which we are 
interested in, and on the observational apparatus which we rely on. The smaller 
scale we pay attention to and the more variety of apparatus we utilize, the closer 
the model Universe constructed from the data approaches the real Universe. 

The author thanks the Ministry of Education, the Government of Japan for 
financial support. He also thanks Inamori Foundation, Japan for encouragement 
as well as financial support. 
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A The zero mode 



In our theory, the zero-mode /o of the Laplacian A is more important than in a 
usual mathematical context: It is directly related to the spatial volume (Eq.( |A2| ) 
below), so that it is a dynamical object just as other modes. Hence it is essential 
to note the basic facts on the zero-mode here for the development of our theory. 

Since A/o = 0, it follows that = / /qA/q^ = — / {dfof ^, implying that 
dah = 0. Thus, 

The zero-mode of A is a constant function, hence there is no degeneracy. 

(Al) 

Now, from Eq.(|l|) for m = n = 0, it follows that 1 = (/o,/o) = IqV on 
account of ( |A1| ). Here V is the {D — l)-volume of the space. ThusQ 



/o = l/VV . (A2) 
Hence J fo^ = VV. From Eq.([^) for n = with Eq.(^), we get 

fmy/=VVSmO . (A3) 



Eq. ([A^) can be represented as /o = (/ y/)^^^"^- Taking the variation of the 
both sides of this equation, we easily get 6fo = —-^^wj2 1 9 " ^9y/- With the help 
of Eq. (^) , we thus get 

5h = -\{g-5g)J, . (A4) 

This formula is of the fundamental importance to develop the perturbation theory 
suitable for our purpose (see Appendix^ and §§ p.3|) . 

On the other hand, taking the variation of Eq. ( |A2| ) directly, we obtain 



5h = -\v-^/^5V . (A5) 

Comparing Eqs. (|A4|) and (^) to each other, we get Eq.(|T^ in 

5V 



B Basic results of the perturbation theory 

The perturbation theory is helpful to to analyze the dynamical evolution of the 
spectra {A„}„=o,i,2,---- present case, however, we need to pay special at- 

tentions to the zero-mode as is explained in Appendix 0. We here derive basic 

® Here we choose the positive squre-root. There is no essential difference even when we 
choose the negative squre-root instead. For instance, Eq.(p^ in §0 remains same. 
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formulas of the perturbation theory, taking care of the zero- mode. (Up to Eg. (P7|) 
below, we mostly follow the argument in Ref.||ll||.) 

We consider an Hermitian operator^ A parameterized by a which is analytic 
in a in the neighborhood of a = 0: 

A = Ao + aAi + c?^2 + ■ ■ • 

=: Ao + a(5„A . (Bl) 

The operator A can be arbitrary although we choose it to be the Laplacian in 
the application. 

Let {(A^j"), |?T.)„)}n=o,i,2,- and {(A^°), |n(°)))}„=o,i,2,- be the set of spectra and 
eigenvectors for A and Aq, respectively: 

Aol^W) = -Ai°V^°^) • (B2) 

We understand that {|n'^°^)}„=o,i,2,--- is the normalized set of eigenvectors. On the 
other hand, we do not specify the normalization of {|^)^}n=o,i,2,-- at this stage 
(see after Eg. (|B7|) below). 

We assume that both A^"^ and \n) ^ are analytic in a in the neighborhood of 
a = 0: 

Ai") = Ar + aA«+a^A(2) + .-- 

=: \f^ab^\n , 
|n)„ = |nW) + «|n«) + «>(2)) + --- . (B3) 

We introduce a projection operator P„ which projects any vector to the sector 
perpendicular to \n'^^^): 

P„:=1_|^(0))(^(0)|^^|^(0)^^^(0)| _ 



With the help of Egs. ( piP and (p3| ), the first eguation in ( P2|) can be repre- 
sented as 

(Af + Ao)|n)„ = -a(5,A„ + 5«A)|n)„ . (B5) 

The L.H.S. vanishes when multiplied by (n^^-'l, indicating that the R.H.S. of 
Eg.(|B5|) is perpendicular to \n^^^). Thus, Eg.(p5D can be represented as 

(AW + Ao)|n)^ = -aP„(5,A„ + 5.A)|n)^ (B6) 

Due to the presence of P„, which removes the zero-mode of the operator A^"-* + Aq, 
Eg.(|B^) can be expressed as 

|n)„ = C„(a)|nW) - P„(^.A„ + 5„A)|n)„ . (B7) 

An + Aq 



It can also be regarded as an Hermitian-operator-valued function of a. 



From here, we need to proceed in a different way from the standard perturbation 
theory in quantum mechanics. Noting that {|?^)Q,}n=o,i,2,--- forms an orthogonal 
set of bases due to Eg. ([B2D , we can interpret Cn{a) = {n^^^\n)^ as the factor 
characterizing the normahzation of {|n)^}„=o,i,2,- - It should be analytic in a 
around a = with C{a) — > 1 as a — > 0: 

= l + «c« + «'ci2) + ... . (B8) 

In the standard perturbation theory in quantum mechanics, the inner-product of 
the states is not perturbed. Thus any normalized state \A) and its perturbation 
S\A) are always perpendicular to each other, so that usually we can set Cnia) = 
1 |]TT|. In our case, on the other hand, the inner-product itself is perturbed 
through the integral-measure when we regard (/„, /„) in Eq.(|I|) as the inner- 
product of the "states" fm and /„. Thus we should tread C„(q;) with care (see 
the argument after Eq.(H)). 

Since (ra(°)|(5c,A„ + 5«A)|n)^ = 0, we get 

Cn{a)6^Xn = -{n^°^KA\n)^ . (B9) 

Eg . (|B9| ) indicates that the both sides should match as analytic functions around 
a = 0. Hence, taking into account Eqs.(|BlD, ( P3| ) and (P8|), we get 

A« = -(n(°)|Ai|nW) , 

A(f) = -(n(°)|Ai|n«)-(nW|A2|n(°))-c«A« , 

(BIO) 

A« = -(nW|Ai|n('-i))-(n(°)|A2|n('-2)) (nW|Ai|n(°)) 

_ Ai-2)\{2) 



Now, the R.H.S. of the second equation in ( |B3[) and the R.H.S. of Eq. 
should match as vector-valued analytic functions around a = 0. Hence, taking 
into account Eq.(|B^), we get 

k 

= , (Bll) 

k 

where 

„(l) ._ J "^^^ mi 9^ 



1^) for m 

' ^(A2)_ + A«/i« +Efc(Ai)^,/iS) for m^n 



for 171 = 12 
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We note that Eg. (|E^) can be expressed as 



Thus 8a\n should be independent of the normahzation factor C„(a). For example, 
with the help of Eg.( pil|) , the formula for A^^^ in Eg. (piO|) turns out to be 

E i^'i-) nki^'i-) kn /A \ (TIT A\ 

—Jo) TioT ~ , (B14) 

k^n ^fc ~ 

which is independent of c^'-* (/ = 1, 2, ■ ■ ■). 

Mostly a is identified with a time-function t in this paper. When we interpret 
the variation SQ of any guantity Q to be a derivative of a corresponding function 
Q{a) at a = 0, only the terms of order 0(a) (guantities with superscript indices 
(0) and (1)) in the above formulas are important. 
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